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io where i and j are item numbers attached to plant values. N is the number of samples of the residual used 
for calculation of the correlation function, and <f>ij is correlation function with respect to the item number j of 
the item number i, which represents auto-correlation function when i=j and cross-correlation function when 

i*j. 

By using the correlation function <t> ti obtained from the above-mentioned equation (3), whiteness 
J5 qualification index ALy defined as the following equation (4) is obtained. 
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l+<j>. 


Then, by making use of the characteristic that ALy depends upon the normal distribution N(0, 1) (having the 
average of 0 and the dispersion of 1), threshold level is determined in accordance with the erroneous 
diagnostic level. 

25 There was conducted a . qualification, such that when the following equation (5) holds, correlation 
:u between the residuals *j and coexists and this has not whiteness. ' ; 

V jALijj > ARL (5) ; : V"" [ ]' V y -[ . ' ' - 

30 In accordance with the above-mentioned diagnostic system, typical abnormality or failure in the plants 
can be identified, but since information i of time constant related to dynamic characteristic of the plant is nor 
included in the above equation (3) for calculating the correlation function fg, identification of abnormality or 
failure in consideration of time constant of the plant to be diagnosed such that influence of abnormality at a 
certain portion would be spreaded to another portion with a certain time delay could not be conducted. 


SUMMARY OF THE INVENTION 


It is an object of the present invention to further improve the plant diagnostic system previously 
proposed and to provide a plant diagnostic system having higher reliability which can conduct identification 

40 of abnormality or failure in consideration of time constant. 

To achieve this object, the present invention is characterized by first means for determining, using the 
coefficient matrix of the auto regressive model, a time constant at which a change in the value of one plant 
variable has the graetest influence on the value of the plant variable and other plant variables; second 
means for calculating a time series of residual differences between an estimated value calculated by using 

45 said plant diagnostic model and an actual plant value obtained from said plant system so as to obtain the 
time series of residual differences; third means for qualifying whiteness of the time series of residual 
differences from said second means and the time constant from said first means so as to produce an index; 
fourth means for determining an abnormality from the index produced by said third means; fifth means for 
determining a failure source from the index and the result obtained by said fourth means; and sixth means 

so for outputting the abnormality and the failure source according to the results obtained by said fourth means 
and said fifth means 

An apparatus according to the invention comprising: 
a plant diagnostic system for a plant having a pluralities of variables, said diagnostic system including; 
means to estimate values of the plant variables using an auto regressive model of the plant to be 

55 diagnosed; means to calculate residuals, these being the differences between the actual values of the plant 
variables and the values estimated using the auto regressive model, means to calculate a correlation 
function on the basis of the residuals, means to diagnose, on the basis of the calculated correlation function, 
whether the plant is operating normally or not, is characterized in that the plant diagnostic system further 
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value is calculated and this value is fed back to the time series of difference calculation means 211. 

The time series of difference calculation means 211 calculates the time series of difference * t (i 
represents an item number of the plant diagnostic model 20) i.e., a time series of difference between an 
estimated value and an actual value in respect of each plant value on the basis of the estimated value from 
5 the plant diagnostic model 20 and a value read from the object 19 to be diagnosed. After the deviation is 
calculated, the time series of difference calculation means 211 calculates an average per each item to 
calculate a residual by subtracting the average from respective deviations to output the time series of 
difference to the time series of difference whiteness qualification means 212. 

The time constant determination means 213 inputs AR model coefficient matrices A(1), A(2), .... A(M) of 
10 the AR model expressed by the equation (1) from the plant diagnostic model 20 to calculate an impulse 
response \\(r) given by the change of a plant value (the item number i) with respect to the change of a plant 
value (the item number j) in accordance with the following equations (6.1), (6.2), (6.3) and (6.4). 

The time constant determination means 213 inputs AR model coefficient matrices A(1), A(2) A(M) of 

the AR model expressed by the equation (1) from the plant diagnostic model 20 to calculate an impulse 
75 response a^r) given by the change of a certain plant value (item number i) with respect to the change of a 
certain plant value (item number j) to obtain sampling time which has the strongest influence on such a 
change. 

For instance, the relationship expressed below holds in connection with the model of M = 2 by the 
equation (1). 

20 

X(s) = ACL) & (s-1) + A(2) £ (s-2) (6.1) 


25<v To what degree the change of the plant value before one sample affects is calculated by the following 
-3 equation. f . • .£V ; - :; ■ ' " ' '; r : * 

% :■ ' . ..... • ' v 

* ... X(s) * A(l) # (s-1) v>.Vv. (6.2) . 

A value of the item number j in the 2K(s) when only a value of.the.jtem number i is set to 1 and values of 
other item numbers are set to 0 in the ^C(s-1) of the equation (6.2)* is an impulse response given by the v 
change of the plant value (item number i) before one sample with respect to the change of a certain value 
35 (the item number j). 

For instance, when the model of M = 2 is employed, i.e., the number of items is 2, the equation (6.1) is 
expressed as follows. 


40 



(6.1) 


The impulse response given by the item number 2 before one sample to the item number 1 (xi(s)) is 
calculated as follows. 

(1) An impulse response ai2(1) given by the item number 2 (x2(s-1)) before one sample with respect to 
so item number 1 (xi(s)) is equal to a value of x 2 (s) in the equation (6.1) when xi(s-1) = 1 and x 2 (s-1) = 0 
in the following equation: 
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A 21 (l) A 22 (l)/ \x 2 (s-l) j 
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20 


25 


given by the change in the plant value (the item number i) before t samples with respect to each plant 
value (the item number j) is expressed as follows. 


T-l 


I. (x) = I A(4) I.(T-A) + x 
1 J^l 1 / A,. (T) 


A li (T) . (6.5) 

2i' 


' A ni< T > 


15 where i, j = 1,2,3 n and r = 1,2 M. and when r = 1, 


I. (T) = I. (1) = / A li (1) 
' A 2i (1 > 


Aji (D 


In addition, time constant ry of the plant value (item number j) with respect to the plant value (item 
« number i) is calculated on the basis of the following definition., 

* . -• •< ./: - : -.v ■> ' ■ ' i. • . ■•• 

30 {rjjautr^ja^r')!; r^r,(r = 1 ,2,...%M),(r' = 1 ^....M)} ^ - ^ . * 

Namely, rij is defined as a value of V. which maximizes the absolute value of a^l), a {j (2), .... a }j (T), .... a fj (M). '<< 
It is considered that, when the plant value of the item number i varies, the time constant thus 
obtained indicates a sampling time duration which has the greatest influence of such a change on the plant 
35 value of the item number j. 

For example, it is assumed that a y calculated from the above equations (6.1) to (6.5) is evaluated as 
follows. 

aij (1) = -3.75 
a,j(2) = -4.71 
40 ajj(3) = 4.08 
aii (4) = 6.10 
ai](5) = 1.24 
a,j(6) = 5.02 
aij(7) = 3.09 
45 au(8) = -2.89 

Because the absolute value of aij(4) is the greatest, the time constant is obtained as rj = 4. 

The time constant r q thus obtained is output from the time series of difference whiteness qualification 
means 212. 

The time series of difference whiteness qualification means 212 receives the time series of difference a 
so output from the time series of difference calculation means 21 1 and the time constant output from the 
time constant determination means 213 to calculate an index for qualifying whiteness of the time series of 
difference. For this, correlation function expressed by the following equation (8) is first calculated. 
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In the thermal power plant thus configured, the above-mentioned three manipulated variables are water- 
fuel ratio (item number 1) of fuel FR, reheater gas damper GD (item number 2) and super heater spray SP 
(item number 3), and the above-mentioned four state variables are addition demand value MWD (item 
number 4), SH outlet temperature MST (item number 5), RH outlet temperature RHT (item number 6) and 
5 main steam pressure MSP (item number 7). 

When the above-mentioned plant variables are input to the plant diagnostic section 21, the time series 
of difference calculation means 21 1 outputs these plant variables to the plant diagnostic model 20 and takes 
therein estimated values of plant variables obtained based on calculation in the plant diagnostic model 20 in 
accordance with the previous outputs. Thus, the time series of difference calculation means 21 1 compares 
io these estimated values with plant values which are input this time from the object 19 to be diagnosed to 
calculate time series of difference ej of each plant variable to output the time series of difference a thus 
computed to the time series of difference whiteness qualification means 212. 

In addition, the time constant determination means 213 inputs AR coefficient matrices A(1), A(2). .... A- 
(M) from the plant diagnostic model 20 to calculate time constant r r j in accordance with the above- 
75 mentioned equations (6.5) and (7) to output the time constant r n thus calculated to the whiteness 
qualification means 21 2. 

The time series of difference whiteness qualification means 212 takes therein these time series of 
difference t\ and time constant r }j to calculate whiteness qualification index AL {j for judging relationship 
between items and abnormality or normarity to output the whiteness qualification index AL fj thus calculated 
20 to the abnormality judgement means 215 and to the failure source determination means 216. At this time, 
by precluding in advance calculation such as whiteness qualification index AUa from the state variable with 
respect to the manipulated variable etc. in connection with the abnormality which actually would not occur 
. in the actual plant, for instance, "the abnormality of SH outlet temperature of the item number 5 has 
influence on SH spray of the item number 3", wasteful abnormality detection operation can be saved. 
25^ Further, when seven plant values are grouped into manipulated variables and. state, variables and the 
^. whiteness qualification index AL y is calculated, e.g.,. by making: use of; the relationship determined so that i 
& = 1 to 3 and j :=:4 to 7. it is also possible to simply manipulated variables which causes the state variable 

to be abnormal. ; •■ '\\ ■ ; > • 

& When the plant is normal, the time series of difference ej. takes random value which becomes non- 
30..' correlated value. The whiteness qualification index AL }j at this time depends upon the normal distribution (0, 
• 1) (having average of 0 and dispersion of 1). ;; 
Thus. . the threshold leveMor judging whether there is correlation: or not can' be determined 1 from the ' 
table showing the normal distribution so that, for example, ARL = 1 .96 at the erroneous diagnostic 
probability of 5% and ARL = 2.56 at the erroneous diagnostic probability of 1%. Namely, when jAL,jj £ 
35 ARL, it is considered that the time series of difference has no correlation. 

The abnormality judgement means 215 inputs ALi 3 indicative of the auto-correlation function of the 
above-mentioned whiteness qualification index AL n to judge the item number i which satisfies the following 
relationship to be abnormal: 

40 jALjjj > ARL. 

The failure source determination means 216 inputs AL {j (i * j) indicative of cross-correlation function of 
the above-mentioned whiteness qualification index AL,j and the item number i judged to be abnormal by the 
abnormality judgement means 215 to judge the abnormality of the state variable (the item number j) to be 
45 based on the failure of the manipulated variable (the item number i) in connection with the manipulated 
variable (the item number i) and the state variable (the item number j) of which both item numbers i and j 
are abnormal. 

It is now assumed that the whiteness qualification index ALy calculated by the whiteness qualification 
calculation means 212 is expressed by the following Table 1 and the threshold value of the erroneous 
50 diagnostic probability of 5% is set to ARL = 1 .96. 


55 


9 


EP 0 216 356 B1 


APC 12 input and therefore they may be input from the APC 12. 

(3) In the identification of the plant diagnostic model, the method was employed to vary at random load 
demand value which was a specified plant state variable thereby to derive mathematical model of the 
object 19 to be diagnosed. However, it is not required that the load demand value is necessarily used in 

5 any plant and it is of course that various plant state variables can be selected in accordance with the 
characteristic of the object 19 to be diagnosed. 

(4) For identification of the plant diagnostic model, random signal such as M-series signal were used. In 
addition, it is possible to identify the plant diagnostic model with another random signal such as pseudo- 
random number. 

70 (5) When not only the signal varying at random but also a step signal or ramp signal is applied to the 
state variable (the load demand value in the example shown in Fig. 3) which has the greatest influence 
on the object 19 to be diagnosed, a method may be employed to test in advance change of the 
whiteness qualification index AL y based on the plant diagnostic model 20, thereby to independently 
determine values of the threshold level ARL per each state variable. 

75 (6) In the above-mentioned embodiment, the threshold level ARL was fixed. However, when disturbance 
which can be measured e.g., large change in load is applied to the object 19 to be diagnosed, for 
reducing the possibility of the erroneous diagnosis, it is possible to allow the threshold level ARL to be 
variable in accordance with the change in the disturbance which can be measured. 
(7) When the diagnostic result is grasped on the basis of elevation of the abnormal level, it is possible to 

20 inform an operator of the approach to the alarm region for paying an attention to the operator and thus 
this can be utilized for preventing accident. 

As described above, according to the present invention, while effecting the control of a plant, the 
diagnosis of the entirety of the plant including its control status can be carried out. In the implementation of 
the diagnosis, the characteristic of the entire plant is first obtained from the plant diagnostic model in terms 
25 of time constant, and then the diagnosis is actually conducted. Accordingly, this makes it possible to 
identify abnormality or failure source in conformity with the characteristic of the object to be diagnosed. In 
addition, since it is unnecessary to test the plant for obtaining the time constant,- the plant can be effectively 
used as compared to other methods. 

30 Claims ' 

. 1. A plant diagnostic system wherein a plant to be diagnosed is expressed as an auto regressive model to 
calculate a correlation function on the basis of a difference between an estimated value of a plant 
variable calculated using said auto regressive model and the actual value of said plant variable, thus to 
35 diagnose said plant, characterized by 

first means (213) for determining, using the coefficient matrix of the auto regressive model, a time 
constant at which a change in the value of one plant variable has the graetest influence on the value of 
the plant variable and other plant variables; 

second means (211) for calculating a time series of residual differences between an estimated value 
40 calculated by using said plant diagnostic model and an actual plant value obtained from said plant 
system so as to obtain the time series of residual differences; 

third means (212) for qualifying whiteness of the time series of residual differences from said second 
means and the time constant from said first means so as to produce an index; 
fourth means (215) for determining an abnormality from the index produced by said third means; 
45 fifth means (216) for determining a failure source from the index and the result obtained by said fourth 
means; and 

sixth means (217) for outputting the abnormality and the failure source according to the results obtained 
by said fourth means and said fifth means. 

so 2. A plant diagnostic system as set forth in claim 1, wherein said first means calculates impulse response 
ajj(r) indicating to what degree a change of a certain plant value (item number i) affects a change of 
another plant value (item number j) from said auto regressive model when auto regressive coefficient 
matrices 

A(1), A(2) A(M) expressed as 

55 
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gegeben sind, um so die Zeitkonstante r n zu berechnen, die den groBten EinfluB hat. 

3. Vorrichtung, umfassend: 

ein Anlagen-Diagnosesystem (21) fur eine Anlage mit einer Vielzahl von Variablen, wobei das 
5 Diagnosesystem folgendes beinhaltet; 

Mittel zur Schatzung von Werten der Anlagenvariablen unter Verwendung eines autoregressiven 
Modells der zu diagnostizierenden Anlage; 

Mittel (211) zur Berechnung von Resten, wobei diese die Differenzen zwischen den tatsachlichen 
Werten der Anlagenvariablen und den Werten, die unter Verwendung des autoregressiven Modells 
70 geschatzt wurden, sind; 

Mittel (212) zur Berechnung einer Korrelationsfunktion auf der Basis der Reste; 
Mittel (215) zur Diagnose auf der Basis der berechneten Korrelationsfunktion, ob die Anlage normal 
lauft oder nicht, 

dadurch gekennzeichnet, dafl das Anlagen-Diagnosesystem (21) ferner ein Zeitkonstanten-Bestim- 
75 mungsmittel (213) enthalt, das die Koeffizientenmatrix des autoregressiven Modells zur Berechnung 
einer Zeitkonstante verwendet, zu welcher eine Veranderung des Wertes einer Anlagenvariablen den 
groBten EinfluB auf den Wert der Anlagenvariablen und die anderen Anlagenvariablen hat, und 

das Mittel (212) zur Berechnung einer Korrelationsfunktion den berechneten Wert der Zeitkonstante 
wie auch die Werte der Reste verwendet. 


20 


Revendications 


25 


1. Systeme de surveillance d' installation dans lequel une installation a surveiller est exprimee selon un 
modele auto regressif pour calculer une fonction de correlation sur la base de la difference entre une 
valeur estimSe d'une variable d'installation calculee en utilisant ledit modele auto regressif et la valeur 
*4 r^elle de ladite variable d'installation, effectuant ainsi la surveillance de ladite installation, caracte rise 
. par. v ' . . \- ••- \ : " 

un premier moyen (213) pour determiner, en utilisant la matrice des coefficients du modele auto. 
v- s $ .regressif, une constante de temps pour laquelle une modification de la valeur d'une variable d'installa- 
30 '.i- tion possede Influence la plus importante sur la valeur de la variable d ? installation et des autres 
: ' variabies d'installation^ 

un second moyen (211) pour cajculer une s^rie temporelle de differences r^siduelles entre une valeur 
estimSe calculee en utilisant ledit modele de surveillance d'installation et une valeur rdelle d'installation 
obtenue a partir dudit systeme d'installation de fagon a obtenir une seVie temporelle des differences 
35 residuelles; 

un troisieme moyen (212) pour qualifier la blancheur de la seVie temporelle des differences residuelles 
a partir dudit second moyen et la constante de temps a partir dudit premier moyen afin de produire un 
indice; 

un quatrieme moyen (215) pour determiner une anomalie a partir de I'indice produit par ledit troisieme 
40 moyen; 

un cinquieme moyen (216) pour determiner une source de defaillance a partir de I'indice et du resultat 
obtenu par ledit quatrieme moyen; et 

un sixieme moyen (217) pour deiivrer en sortie I'anomalie et la source de defaut selon les resultats 
obtenus par ledit quatrieme moyen et par ledit cinquieme moyen. 


2. 


Systeme de surveillance d'installation selon la revendication 1, dans lequel ledit premier moyen calcule 
la reponse impulsionnelle ay (t), indiquant dans quelle mesure la modification d'une certaine valeur 
d'installation (element numero i) modifie une autre valeur d'installation (element numero j) d'apres ledit 
modele auto regressif lorsque les matrices des coefficients auto regressifs 
so A(1 ), A(2) A(M) exprimes par 

~ M 

X(s) = 2 A(ia)-X(s-m ) 

55 m=l 

sont donnees pour calculer ainsi ladite constante de temps r r) qui possede I'influence la plus 
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